Some exact infrared properties of gluon and ghost propagators 
and long-range force in QCD 
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We derive some exact relations in Landau gauge that follow from a cut-off at the Gribov horizon 
which is then implemented by a local, renormalizable action involving auxiliary bose and fermi 
ghosts. The fermi ghost propagator is more singular than 1/k 2 at k = 0, and the relation ao+2aa = 
(D — 4)/2 holds between the infrared critical exponents of the gluon and ghost propagators D(k) 
and G(k) in D Euclidean dimensions. Finally, in D Euclidean dimensions, there is a long-range 
force, transmitted by the propagator of the auxiliary bose ghost that corresponds to a linearly rising 
potential with tensor coupling to colored quarks that is proportional to the renormalization-group 
invariant g 2 D(k)G 2 (k). A comparison with numerical results is discussed. 

PACS numbers: 12.38.-t, 12.38.Aw, 12.38.Lg, 12.38.Gc 



I. INTRODUCTION 



The present study is motivated by recent numerical in- 
vestigations on very large lattices [l|, [H, H, H, [f| , in which 
it is reported that, in D = 3 and 4 Euclidean dimensions 
in minimal Landau gauge, 1 the gluon propagator D(k) 
approaches a finite value at k = 0, while the ghost prop- 
agator G(k) behaves like 1/k 2 . These properties are in 
contradiction with a scenario developed by Gribov Q and 
elaborated by the author 0, [1| according to which the 
gluon propagator D(k) is suppressed and even vanishing 
at k — 0, while the ghost propagator G(k) is enhanced at 
k = compared to 1/k 2 . On the other hand numerical re- 
sults in D = 2 Euclidean dimensions are consistent with 
this scenario [T3| , and moreover in a detailed numerical 
study of Landau gauge in D = 3 Euclidean dimensions, 
it was found that if one gauge fixes to the extent prac- 
ticable to an absolute — instead of merely a relative — 
minimum of the minimizing functional, then results are 
also not inconsistent with the Gribov scenario It 
was also observed that in 4 dimensions it would be pro- 
hibitively costly to seek an absolute minimum on lattices 
of sufficiently large size . 

In this unclear situation, it seemed worthwhile to see 
what exact properties could be derived by analytic meth- 
ods, if the functional integral is cut-off at the Gribov hori- 
zon. We start by reviewing previous results which are 
scattered in different articles. We then obtain as exact 
infrared properties in D Euclidean dimensions (i) that 
the ghost propagator G(k) is more singular than 1/k 2 , 
(ii) that the equation fl2l. Hij, fl4| . 



and qg are defined by, 
D(k) 



(k 



G(k) 



(k 



2U + QG 



(2) 



and (iii) that there is a long-range linearly rising 
renormalization-group invariant potential that has ten- 
sor coupling to colored quarks. 



II. ABSOLUTE LANDAU GAUGE 

The absolute Landau gauge is defined by gauge trans- 
forming to the absolute minimum of the functional [l5| 



Fa(9) = \\ 9 A\\ 



(3) 



with respect to all gauge transformations g, where 9 A = 
g~ l Ag + g _1 dg is the gauge transform of A, and \\A\\ is 
the Hilbert norm of A. From stationarity at a relative 
or absolute minimum, it follows that A is transverse, d ■ 
A = 0, so the fundamental modular region A, which is 
defined to be the set of absolute minima of the minimizing 
functional .FaQ?) on each gauge orbit, satisfies 



A = {A:d-A = 0,\\A\\ < \\ 9 A\\}, 



(4) 



for all gauge transformations g. 

The set of relative minima, called the Gribov region Q, 
is larger than the set of absolute minima, A C fl, and may 
be characterized as the region, where A is transverse and 
the Faddeev-Popov operator is non-negative, M(A) > 0, 
which means all its eigenvalues are non-negative, 



a D + 2a G = {D- 4)/2, 



(1) 



holds, where the gluon and ghost infrared exponents ao 



1 Here the gauge is fixed at a relative minimum of the minimizing 
functional 151 . 



Q = {A : d ■ A = 0, M(A) > 0} 
The Faddeev-Popov operator is defined by 
M ac (A) = -D^d^ = -d 2 S ac - gf 



(5) 



(6) 



where D°f is the gauge-covariant derivative and a, 6, c are 
color indices, with a = 1, ... N for SU(N). The positivity 
property holds because at a relative or absolute minimum 



2 



the matrix of second derivatives is positive which, for the 
minimizing functional ([3]), is the operator M{A). 

Recall that the fundamental modular region A enjoys 
the following simple properties [l6j . (i) The configuration 
A = is an interior point of A. (ii) A is bounded in 
every direction, which means that if A £ A then for some 
sufficiently large positive number A > 1 the configuration 
XA ^ A. (iii) A is convex, which means that if A\ G A and 
£ A, than for positive numbers a > and j3 > with 
a + [3 = 1, then the linear combination, aA\ + (3A2 G A, 
lies in A. The Gribov region 17 enjoys the same three 
properties [l7| . 



III. CUT-OFF FUNCTIONAL INTEGRAL 

In the absolute Landau gauge the partition function is 
expressed as an integral over A, 



discrepancy noted in the Introduction between the prop- 
erties we shall obtain here and numerical results obtained 
by fixing to a relative minimum. 

The advantage of integrating over the Gribov region f2 
is that it is explicitly described as the set of transverse 
configurations A for which the "horizon function" H (A) 
has the upper bound 

n = {A : d ■ A = 0, H(A) < qV}, (10) 

where Vis the Euclidean volume, and 



q = (N 2 — 1)D, 



(11) 



is the number of components of the gauge field A b ^. 
Points A on the boundary dfl of the Gribov region sat- 
isfy H(A) = qV. The horizon function is given explicitly 
by 0, 



Z A = N J dA 6{d ■ A) exp[-S(A)}, 



(7) 



H = g 2 {A 1 M- 1 A) = J d D x h{x) 



(12) 



Here 



where h{x) is the density 



S{A) = -lndet(M(A))+S YM (A))-lndet(P+m q ), (8) h ( x ) = 9^ f a ° c ^{x) I d D y (M l ) c *{x, y; A) f ade A^y) . 



is an effective action that includes the Faddeev-Popov 
and quark determinants, and Sym is the Euclidean 
Yang-Mills action. The Faddeev-Popov determinant, 
dct(M(A)) : is positive for A in A because A C fi, and 
M(A) is a positive operator in il. 

We do not have an explicit analytic expression for the 
boundary of the fundamental modular region A. How- 
ever we know that it is contained in the Gribov region £1 
and that part of their boundary is common. This part in- 
cludes Abelian configurations and, on the lattice, center- 
vortex configurations which, in some scenarios, are the 
dominant configurations [lg]. It has been argued that 
for purposes of the functional integral, the fundamental 
modular region A and the Gribov region Q arc equiv- 
alent [ijj In the present article we take this as a 
working hypothesis, and calculate Euclidean averages by 
integrating over instead of A, with partition function 



Z n = N dA 6(d ■ A) cxp[-5(A)] 



(9) 



We shall deduce some exact properties of the gluon and 
ghost propagators, and discuss the validity of the results 
obtained in the concluding section. 

When comparing with numerical determinations of the 
gluon and ghost propagators, it should be kept in mind 
that the lattice analog of integrating over f2 is not the 
same as numerical gauge fixing to a relative minimum 
by a particular algorithm. For although such an algo- 
rithm does give a gauge with a probability distribution 
that lies within f2, this gauge is, in principle, algorithm- 
dependent. This distinction may be the origin of the 



(13) 

This expression for H(A) resembles the formula for the 
T-matrix in potential scattering at zero energy, T = 
V — VH~ 1 V, where H is the Hamiltonian. It was ob- 
tained by a quantum mechanical-type calculation that 
consists in summing the perturbation scries for the low- 
est non-trivial eigenvalue of M(A). Note that H(A) > 
is positive for A in f2. The cut-off is non-perturbative be- 
cause pcrturbatively we have H(A) — 0(g 2 ), while the 
cut-off occurs at H{A) = 0(1). 

We write the partition function Zq as 

Z n = N J dA 5{d ■ A) 6(qV - H) exp[-S{A)], (14) 

where the restriction to is effected by the ^-function. 



IV. BOLTZMANN-TYPE DISTRIBUTION 

In this section we shall convert the functional inte- 
gral Zq with a cut-off at the boundary of f2 into a 
Boltzmann-type distribution. 

For this purpose we represent the ^-function by its 
fourier decomposition, 



8(qV -H) = J 



duj exp[iuj(qV — H)] 



2iri oj — ie 

and perform the A-integration first, so 



Z n = — 



diu Z(iu>) 



2ir e + ioj 



(15) 



(16) 
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Here 

Zfr) = N J dA8{d- A) exp[-5 + j(qV - H)], (17) 

is a Boltzmann-type partition function that depends on 
a thermodynamic parameter 7 that is the analog of 
j3 = 1/kT, and that may be complex. The denomina- 
tor in (fT6|) is analytic in the lower half w-plane defined 
by u — > u — vy, or iu — ► 7 + iw, where 7 > is posi- 
tive. The partition function Z(j + iw) is also analytic in 
z = j + iuj for 7 > 0. Indeed cxp[iu(qV — H)] is an entire 
analytic function of u, and moreover in the lower half 
plane the factor exp[— (7 + iu)H(A)] only improves con- 
vergence for the A- integration because H(A) > 0. Thus 
we may deform the contour of integration in (|16p into the 
lower half-plane, 



Zr 



du exp[jy(7 + iu) - ln(7 + iu)], (18) 



where the free energy is defined by 

W(j + iu) = ln[Z(7 + iu)]. 
It is a bulk quantity, 

W(rr) = Vw{j) + o(i). 



(19) 



(20) 



Here V is of the order of the total number of degees 
of freedom which, with a lattice cut-off, is of the order 
of the volume of the lattice in lattice units, V jaP . We 
are interested in the limit V/a D — > 00. which we write 
as V — > 00. 

For large V, the exponent W(j) — 0(V) in flH]l is 
large, so the contour integration over u may be evaluated 
by the saddle-point method. To leading order in V in the 
exponent, the term ln(7 + iu), which is of order 1, may 
be neglected, and the saddle point. j s , if it exists, is the 
solution of 



dW( ls ) 
(?7 



0. 



(21) 



This saddle-point condition is a gap equation that re- 
lates 7s to the mass scale Aqqd [3]- 

In the neighborhood of the saddle-point, W has the 
expansion 



W(-y a + iu) 

where, by (fTT)) . 

d 2 W( ls ) 



W{ ls ) 



d 2 W{ ls ) uP_ 
<9 7 2 2 



<9 7 2 



= (H 2 ) j3 - (H)l > 0. 



(22) 



(23) 



The subscript refers to an expectation-value calculated 
in the Boltzmann-type ensemble (jTTJ) namely, 

(X) 7 = Z~\i) (24) 

x N J dAS(d- A) exp[-S + ^{qV - H)] X. 



The quantity on the right of (|23|) is the variance of the 
horizon function H. It is positive and of order V, so the 
width of the peak in the w-integration is of order V^ 1 ^ 2 . 
Thus the saddle-point approximation becomes exact in 
the limit V — > 00, and we obtain the basic relation 



Zq = Z( ls ). 



(25) 



This equality represents symbolically the equality of ex- 
pectation values calculated with the cut-off distribution 
Zq and the Boltzmann-type distribution Z(p( s ) as in (|24p . 

A. Remarks 

(i) We would have gotten the same result if the 9- 
function, 9(qV — H) in (fl4|) . were replaced by the 5- 
function 5(qV — H). Indeed we have seen that the term 
\n(^ + iu) in (|18p may be neglected to leading order in V. 
This term comes from the denominator u — ie in (|15| , so 
we have in fact replaced the denominator by 1. This is 
the same as replacing 9(qV — H) by the <S-function, 



8(qV-H) = J ^expMfV-H)], 



(26) 



whose fourier decomposition has no denominator, 
(ii) Since 8(qV — H) is entirely concentrated on the 
boundary H = qV of the fundamental modular region, 
it follows that at large V the probability gets concen- 
trated on the boundary, (iii) The microcanonical en- 
semble 5{qV — H), and the equivalent cut-off ensem- 
ble 9{qV — H), are also equivalent to a canonical or 
Boltzmann-type ensemble (|17[) . where the parameter 7 

For, by (fl"7|) . the saddle-point 

= 0, with W(^) = lnZ(7), is equivalent 



is the analog of (3 = 2 



condition, ^Ili — 



to (qV - H) 1 = 0, or 



(H) 



qV, 



(27) 



where the subscript is defined in (|24| . Thus the "hori- 
zon condition" H = qV, which would be imposed by 
6(gV — H) is satisfied in the mean in the Boltzmann- 
type distribution. Moreover, since H(A) is a bulk quan- 
tity, its fluctuations are of relative order l/V 1 / 2 which is 
negligible. Thus the saddle-point condition is the horizon 
condition. 



V. LOCAL RENORMALLIZABLE ACTION 

We wish to represent the partition function Z(j), 
which contains the non-local action S + 7 (if — fV) 



2 Here we use the terms "microcanonical" and "canonical" by anal- 
ogy with the corresponding distributions in statistical mechanics 
to which they are mathematically, but not physically identical. 
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as a functional integral with a local action. We use 
the Faddeev-Popov identity as usual to rewrite S(d ■ 
A)detM(A) as an integral over additional fields, with 
the local Faddeev-Popov action 



dAS(d-A) det M(A) exp(-Sy M ) 

= J dAdbdcdc exp(-SW), (28) 



where Sfp = J d D x Cfp is the local Faddev-Popov ac- 
tion in Landau gauge 



pacd 



Cpp — — F, 



idubAu - ducDuC, 



and 



F 



(29) 



(30) 



where (A^ x A v ) a = f abc '-A h A c v . The Lagrange-multiplicr 
field b imposes the Landau gauge condition d ■ A = 0. 
Without loss of generality we ignore the quark action 
which plays no role in the discussion. 

As a second step, we similarly eliminate the non-local 
term j(H — qV) in the action (| LT[) by an integral over 
auxiliary fields and a local action. This is accomplished 
by the identity @,H, 



cxp[— j(H — qV)] = / d(j)d<j)dujdtjj exp(—S a ux — S 7 ), 



that is easily verified by integrating out the auxiliary 
fields by Gaussian integration. They consist of a quar- 



tet of auxililary bose and fermi ghosts, 



and 



and corresponding anti-ghosts and ui ab that carry a 
Lorentz index and a pair of color indices. The term on 
the left, jH = jg 2 (A, M~ 1 A), is obtained by completing 
the square in the exponent in the Gaussian integration 
of the bose ghosts <fi and <p. This produces an unwanted 
l/det q (M) in the denominator which is then cancelled 
by det 9 (M) in the numerator that comes from the inte- 
grating out the fermi ghosts to and co. The total local 
action is given by 



S = 



c 



where 



and 



£ — £fP "t - £aux 



(32) 



(33) 



(34) 



The gauge-covariant derivative and the Lie commutator 
act on the first color index only, while the second color 
index is mute, thus {D X (f> l ,) ab = d A < b + g{A x x ^) afc 



where {A\ x cf> ti ) ab = f 
the thermodynamic parameter 7, 



A c x 4>^. The last term involves 



= 7 1/2 



D\(4>\ - 4>\) - g{D\c x u\) 



79 . (35) 



If the term £ 7 were absent, the integral over the bose 
and fermi ghosts would produce cancelling factors of 
det M, and the action would be strictly equivalent to the 
Faddeev-Popov action. 

All terms except £ 7 are of dimension 4 and respect 
a BRST symmetry transformation that acts on the 
Faddeev-Popov fields according to 



sA 



sc. 



ib; 



sc = —(g/2)c x c 

sb = 0, (36) 



and acts trivially on the auxiliary ghosts, 



s 

SCO 



ab 
I' 
ab 
I' 



ab . 

H ! 

ab . 
I' ' 



SCO 

s 



ab 
ab 



: 

0. 



It is nil-potent, s 2 = 0. The term 

Caux = sd X UJ ab {D x ^) ab 

= d x r b {D^r b 



(37) 



(38) 



d x Qf[ (D x uj,r b + (gD x cx^) 



a b 



is s-exact. 

This symmetry is softly broken by the term in £ 7 that 
is of dimension 2, so renormalizability is preserved 0, 



(31) |2l|, [22j. [The terms in C aU x and £ 7 that contain the 



Faddeev-Popov ghost c and the auxiliary ghost Co play 
no dynamical role because there is no corresponding term 
that involves c and to, so these terms may be eliminated 
by a suitable translation of to.] 

The local, renormalizable action S(j) is an "extended" 
action in the sense that it depends on the arbitrary pa- 
rameter 7. It corresponds to a gauge theory only when 
the saddle-point condition (|27|) is satisfied. 



VI. INFRARED LIMIT OF FERMI GHOST 
PROPAGATOR 

In this section we show that the horizon condition (|2"T|) 
is equivalent to the condition that the fermi-ghost prop- 
agator G(k) be more singular than 1/k 2 at k = 0, 



Hm[rG(fc)]- 1 = 0. 



(39) 



Here 



S ab G(x -y) = (c a (x)c»(y)) = ((M ) ab (x, y; A)) (40) 

is the fermi-ghost propagator, and G{k) is its fourier 
transform. This form of the horizon condition makes 
manifest that it is multiplicatively renormalizable. This 
result was obtained previously 0, 0|. We give here a 
simpler, more direct proof. 
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As a first step we develop a useful expression for the 
horizon or saddle-point condition ([27)) in terms of local 
fields. By translation invariance and (|12p . the horizon 
condition may be written 

(h(0))=q, (41) 

where the expectation value is calculated in the ensem- 
ble defined by the action S. The kernel of the inverse 
Faddeev-Popov operator is the expectation value of the 
fermi-ghost correlator at fixed A, 

(c a (x)c b (y)) A = (M- 1 r b (x,y;A). (42) 

This, and expression (fT3| for h(x), yields 

(h(0))=g 2 ( r bd Al(0)c d (0) J d D y r*A*(y)c?(y) ). 

(43) 

In terms of the 4-point function 

A**?(x,z,y,w) = (Al(x) c d (z) Al{y) c») (44) 
this reads 

(MO)) = 9 2 f abd f ace J d D y A^ e (0, 0, y, y). (45) 

The integral over y projects onto the zero-momentum 
component, and we have 

(HO)) = 9 2 f abd r ce A^f(fc)U=o. (46) 

where 

A bd l e (k) ee J d D x exp(*fc • y) A^f (0, 0, y, y) (47) 

is a fourier transform. The 4-point function (|44f has the 
cluster decomposition into connected parts, 

A bd ^{x,z,y,w) = (A^x) Al(y))(c d {z) c») 

+A con bd l e (x lZl y,w), (48) 

which gives 

A^ e (fc) = [D% G de )(k) + A con bd ^(k), (49) 

where 

D%{x-y) = (Al{x)Al{y)) (50) 
is the gluon propagator. Thus we obtain finally 

(Mo)) = g 2 f abd f ace [{D% G de ) + A con bdc ; ](k)\ k=0 , 

(51) 

and the horizon condition q = (h(0)) is expressed by 

q-g 2 f abd f ace [{D% G de ) + A con bdc ; ](k)\ k=0 = 0. (52) 

We now compare this form of the horizon condition 
with the Dyson-Schwingcr equation (DSE) for the ghost 



propagator G, which we write symbolically to avoid a 
profusion of indices, 

G _1 (Jfe) = k 2 ~gfk ^(loop)(fc) 

= k 2 -gfk Y, G cc D A^ c ^ rc (k). (53) 

This expression is summed over loops namely over 
fields $ that have non-zero propagators with A namely 
$ = A, <j> or <f). Here k 2 is the tree-level term, gfk 
is the tree level vertex, where / represents f abc which 
is contracted as in (|5ip . The external ghost momen- 
tum k factors out of the loop integral because of the 
transvcrsality of the gluon propagator in Landau gauge, 
{k - p) li ,D flv (p) = k^D^ip). 

The DSE for the 3-point vertex reads 

r c ,*, e = kfg + (loop) kfg (54) 

where the external anti-ghost momentum k again factors 
out of the loop. When this is substituted into the ghost 
DSE, both the ghost and anti-ghost momenta factor out, 
and it reads 

G-^k) = k^(k)k„ (55) 

where 

iWfc) = 6^ - gf [(D^G) + A' c „]{k) fg. (56) 

Upon examining the sum of the loop terms, one finds 
that A' contlv (k) is the sum of all terms in the skeleton ex- 
pansion of the connected 4-point function defined above, 

A con ^ v {k) = A' contlv (k) + A S i ng4lu (k), (57) 

except for the piece that is connected by a single ghost 
line G(k), 

A smg (k) ee (GDT)^(k) G(k) (TGD)u(k), (58) 

which is missing. Thus, by comparison with (|5ip we ob- 
tain the interesting relation 

q-(h(0)) = (n + A sing )™(k)\ k=0 . (59) 

The term A S i„ g ^° (k)\ k= o is somewhat ambiguous, be- 
ing of the form 0/0, but we will determine its value at 
k = by Lorentz invariance. By Lorentz invariance 
(GDT) fi (k) = h(k 2 )k^, and the vector fc M has no direc- 
tion at k = 0, so each of the factors in (|58f vanishes at 
k = 0. Thus the missing term in fact vanishes at k = 0, 
and we obtain 

g-(Mo)) = n-(fc)| fe=0 . (60) 

Thus the horizon condition reads 

IC(fc)U=o = 0. (61) 
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This form of the horizon condition is free of the singular 
term A s j„ s (0). 

By Lorentz and color invariance we have 

n^(fc) = [ a(k 2 )S^ + b(k 2 )k„k v } S cd . (62) 

Here a(k 2 ) and b(k 2 ) arc regular at k — if n M „(fc) is reg- 
ular, as appears from its skeleton expansion. At small k 
the second term is negligible compared to the first. Thus 
at small k we have 

U flv (k) w aS^, (63) 

so at small k 

MW*) fc »- » k 2 Tl^(k)/D, (64) 

which gives 

lim i^Gik)]- 1 = n w (0)/D = 0, (65) 

k— >0 

by (1551 and (jrJTj) . which establishes (|39|) . as asserted. 

VII. LONG RANGE OF BOSE GHOST 
PROPAGATOR 

We calculate the bose-ghost propagator in the infrared. 
For this purpose we diagonalize the action by writing the 
bose ghosts as 

<j)= (U + iV)/V2; 4>=(U-iV)/V2, (66) 

where U and V are real fields. The part of the action 
that contains <f> and cf> decomposes into a sum of terms 
that involve U and V separately, 

£-4>,4> = &v + Ay + Am (67) 

where 

£v = \dxV^ h {D x V^) a \y) + iy/2^/\D x V x ) aa (68) 

Ay = \dxUf(D x U,r b (y). (69) 
The mixed term 

C m = {zg/2)d x {r h X d K l ) A l ( 7 °) 

vanishes by partial integration when A is transverse, 
which makes M(A) hermitian. This mixed term is elim- 
inated in the local formalism by a shift of the Lagrange 
multiplier field, 

b = b'- {g/2)f abc V bd U; d . (71) 

The [/-field does not mix with A, and the U-U 
propagator only appears in closed ghost loops, just 
like the fermi ghosts. However because of the term 



iV2-f 1/2 f abc A b x V^ a that is contained in the last term of 
Cv, the field V x does mix with the gluon field A x in the 
sense that the mixed propagator Dav, which is given in 
the Appendix, is non-zero. 

We now calculate the long-range part of the V- 
propagator 

D vv ^f{x -y) = «(*)y; d (y)>. (72) 

The V^-ficld appears at most quadratically in Cv, and we 
obtain a formula for the ^-propagator, by integrating the 
F-field out of the action Sv = / dPx Cv- This is done by 
Gaussian integration after completing the square by shift- 
ing the integration variable V x b by iy /r 2g-f 1/2 (M~ 1 A\) ab , 
where 

(M- 1 A x ) ab (x) = J d D z(M- 1 r c (x,z;A)r db A d x (z), 

(73) 

with the result 
D V vf c l f^-y) = Sx fl S bd 6 ac G(x-y)+Fff(x-y), (74) 

where 5 ac G(x—y) = {{M~ l ) ac (x, y, A)) is the fermi-ghost 
propagator, and 

V) = -27.9 2 ( (M- l A x ) ab (x) (M- l A,f d {y) ). 

(75) 

We now show that there is a very long range compo- 
nent to the F-term in the V propagator. By (|42|) we may 
replace M _1 by a fermi-ghost pair, so 

F^ d {x-y) = -2 lg 2 J d D zd D wKl b l d {x lZl y,w) (76) 

where 

Kf; d (x,z,y,w) ee ( C »c e (z) r b A{(z) (77) 
xuf (y)u)f (w) f hld A^(w) ) 

is a 6-point correlator. Here we have arbitrarily cho- 
sen one of the equal components of the lo-uj propagator 
which are each equal to the propagator of the Faddeev- 
Popov ghost (uif 1 (y)Q bl (w)) = (c a (y)c b (w)) . It may not 
seem like progress to express the 2-point function (VV) 
in terms of the 6-point function K, but in fact the cluster 
decomposition of K enables us to determine its leading 
infrared behavior. The leading term in its cluster decom- 
position is the product of three connected pieces, 

K 3 ff(x,z,y,w) = (c a (x)c E (z)) f^ b {A^A^w)) 
xf hld (u;f(y)^\w)) 
= G(x-z)f aeb D x ^z~w) 

xr d G(w-y). (78) 

The integral over z and w converts this into a convolution 
which, in momentum space, is the product of propaga- 
tors, 

F 3 ,x b f(k) = -2 7ff 2 G{k)r b D x ^{k)r d G{k). (79) 
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This is highly singular at small k. For example if G 
and D were canonical, free, masslcss propagators (which 
they are not) the product would be l/(fc 2 ) 3 . This is the 
most infrared singular part of the ^-ghost propagator, 



D 



w,A?(ft) * -2ig 2 G 2 (k)D Xtl (k)rT d , 



(80) 



at small k. 

We now examine the properties of this asymptotic in- 
frared propagator. The gluon propagator 



D Xll {k) = D(k)P^(k) 



(81) 



is transverse in Landau gauge, where P^ v {k) = 5^ — 
k^ky/k 2 is the transverse projector. We write the long- 
range part of the bose ghost propagator as 

Dy^ffik) « Px^k)r b r d DvvAk), (82) 

which contains transverse and adjoint projectors. The 
long-range scalar propagator is given by 



D VVilr = -2 19 2 G 2 {k)D{k). 



(83) 



The combination g 2 G 2 D is a renormalization-group in- 
variant in Landau gauge [IH, G3 • This is a consequence 
of the factorization of the external ghost momenta in the 
ghost-ghost-gluon vertex, from which it follows that this 
vertex is finite and does not require renormalization in 
Landau gauge, so 



Z 1 ' 2 7~ - 

J 9^3 



(84) 



The combination g G (k)D(k) renormalizes by 
Z 2 Z^Z 2 = 1. It has been proposed that this com- 
bination be used to define a renormalization-group 
invariant running-coupling in Landau gauge fl2| , 



(k) =g 2 (k 2 yG 2 (k)D(k)/(2ir). 



(85) 



The infrared behavior of the gluon and ghost propaga- 
tors is described by infrared exponents ([2]). In terms of 
these we have 



DvvAk) ~G 2 {k)D(k) 



1 



(86) 



0.2)3+aD+2a G ' 

which gives for the infrared power law of the bose ghost, 

1 



D 



VV,lr 



(k) 



( fc 2)l + D/2' 



(87) 



in D Euclidean dimensions, as follows from |T]). 

[We sketch the derivation of equation (fTJ). It is the 
same as in Faddeev- Popov theory, when the horizon con- 
dition [fc 2 G(/c)] _1 |k=o = holds [H, [ll] , because the 
DSE of the fermi ghost is the same. The horizon condi- 
tion iLj„(0) = and the fermi-ghost DSE (|55"j) together 
imply that the ghost propagator satisfies the equation 



G~ 1 (k) = fc AI [n^(fc) - n M „(o)]fc„, 



(88) 



or, by flU, 

G-\k) = k^gf[(D^G)(Q) - (D^G)(k) 



li/u/ (0) - A' con ^ u (k)]fgK. (89) 



The terms involving D^ U G represent a loop integral 
which, for small values of the external momentum k, 
may be evaluated using the asymptotic infrared expres- 
sions (|5J), 



G-\k) 



d D 



kl- 



ip 



2\l+a L 



1 



1 



[{p + kf 



)211+a G r(p)211+a G 

(90) 

On the left-hand side we have G -1 (fc) ~ (k 2 ) 1+aG , 
whereas on the right-hand side p scales like k, so the 
right-hand side is proportional to (fc 2 )" 1 ~ [ * D ^ aG+D / 2 . 
Upon equating like powers of momentum on both sides 
we obtain ([1]). The correction terms involving A' con do 
not change this power counting [l4| ■] 

We contrast the exact power law ([57)1 for Dyv.ir with 
the semi-perturbative expansion that starts from the 
quadratic terms in the local action 5*. (We call it 'semi- 
perturbative' because the non-perturbative horizon con- 
dition is imposed after calculating to a given order in 
perturbation theory.) Indeed the transverse part of the 
U-propagator, calculated from the quadratic parts of S, 
is given by Gribov's expression, 



Dvv = Daa 



k 



(k 2 y 



(91) 



This vanishes like k 2 at k = 0, whereas we the exact 
result ([87]) diverges like l/k 6 in 4 dimensions. Thus semi- 
perturbative theory can be a very bad guide. 

If the power law of Dw.ir, cq (|87|) . described a poten- 
tial, it would be over-confining by l/k 2 . Indeed in the 
non-relativistic limit k — ► k, it behaves like l/k 2+ -° = 
l/|k| 3+s , in S = D — 1 spatial dimensions, whereas a 
linearly rising potential in S spatial dimensions behaves 
like l/|k| 1+s . However we must also consider how the 
bose ghost V couples to quarks. 



VIII. BOSE GHOST AS CARRIER OF 
LONG-RANGE FORCE 

The ghost V does not couple directly to quarks, but 
there is a non-zero q-q-V vertex because the mixed prop- 
agator Dva is non-zero. For example, the triangle dia- 
gram 



g 3 X a ^S qqlu X b Daa^x D AV X e f acd k x 



(92) 



contributes to the q-q-V vertex. Here A a are the Gell- 
Mann matrices, and gf acd k\ is the tree level V-V-A 
vertex. Consequently a F-quantum may be exchanged 
between quarks. 

The external ghost momentum k\ that appears here is 
a general feature of the Landau gauge. It appears at ev- 
ery vertex with an external ghost line, as explained below 
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equation (|53|) . This happens for all ghosts including the 
y-ghost because in the local action £y, the only vertex 
in which the V -field appears is (g/2)d x V° b (A x x V^) ab . 
Thus, when the V propagator Dyv(k) connects two 
quark lines, there is a factor k at each vertex, so the 
effective propagator is of the form 

which rises linearly in position space. However because 
there are two Lorentz indices, the coupling is of tensor 
rather than vector character. 

It would take a non-perturbative, dynamical calcula- 
tion to determine the q-q-V vertex. We do not attempt 
this here. However we can assert that if this vertex is 
finite at k = 0, after factorization of the external V- 
ghost momentum k, then the effective propagator be- 
haves like l/k 1+s . As an example of what the vertex 
could be, we may determine its form at small k if we also 
assume, for simplicity, that there are no powers of quark 
momentum at the vertex. The vertex contains one factor 
of kfj, , as we have seen, and is contracted with the vector 
field V v . The long-range propagator term we have found 
is in the adjoint representation, so this vertex is of the 
tensor form T^k^. Since we assumed that there are no 
further powers of momentum, it has the Dirac structure 
q(A5^ + Ba^)\ a q fc M . The term S^k^ = k v is longitu- 
dinal, and thus orthogonal to the transverse long-range 
propagator we have just found. This gives a q-q- V vertex 
of the form 

q<J^\ a q k„. (94) 

With this vertex, the net power of the exchange of a V 
particle between a pair of quarks is given by 

^A a 5 {k 2^ k + D/ 2 (95) 

In the non-relativistic limit, with spatial dimension S = 
D — 1, this reduces to 

qa^q SijV ^~* h * j q'^XW- (96) 

In position space this is a linearly rising, spin-dependent 
potential. 

IX. CONCLUSION 

We have found that it is possible to derive exact and 
unique infrared properties from Faddeev-Popov theory 
with a cut-off at the Gribov horizon. The uniqueness 
contrasts with the situation of the DS equations of the 
standard Faddeev-Popov action whose solution is not 
unique because one may cut off the functional integral at 
any one of the Gribov horizons without changing the DS 
equations [l3|, and a one-parameter family of solutions 



has been found [23[. It has been proposed to overcome 
this non-uniqueness by imposing an additional condition 
[k 2 G(k)]- 1 \ k=0 = 0. 

One old lesson learned anew here is that it is dangerous 
to make approximations when non-perturbative effects 
are significant. In particular the "free" propagator of the 
bose ghost V® b that results from the quadratic part of 

the action S is given by k 2 /[(fc 2 ) 2 + m 4 ] which vanishes 
like k 2 in the infrared, whereas we have found that its 
exact behavior in the infrared is 1/fc 6 . The Gribov pole 
which occurs at the unphysical value (fc 2 ) 2 + m 4 = 
or k 2 = ±im 2 could turn out to be an artifact of the 
free propagator which, as we have seen, can be highly 
unreliable. The long-range propagator of the bose ghost 
may offer a solution to the problem in Landau gauge of 
the origin of a long range force between quarks. 

It is encouraging that the exact properties derived here 
are consistent with the numerical results of Maas [ll[ that 
are obtained by fixing to the absolute Landau gauge to 
the extent practicable. He has proposed a possible expla- 
nation for why numerical results, described in the Intro- 
duction, are different in minimal and absolute Landau 
gauges. The exact properties we have found disagree 
with the numerical results in the minimal Landau gauge 
in D = 3 or 4 Euclidean dimensions. However, as noted 
in sec. II, the procedure adopted here, which consists of 
integrating over the Gribov region, is not the same as 
a minimal Landau gauge which is algorithm dependent. 
This may explain why the properties we find agree rather 
with the absolute Landau gauge. 

The working hypothesis on which the present article is 
based — that for purposes of the functional integral the 
Gribov region O and fundamental modular region A are 
equivalent — does not have as firm a basis as one might 
wish. Nevertheless the internal consistency of results de- 
rived from this Ansatz is impressive. Remarkably, the 
cut-off at the Gribov horizon which is non-local in the 
gauge field turns out to be implementable by a local 
action that is renormalizable, and moreover the renor- 
malization constants are the same as in Faddeev-Popov 
theory 0). This happens because the modifying term is 
soft, of dimension 2. The Gribov parameter 7 does not 
have an independent subtraction, and the derived condi- 
tion [k 2 G(k)]~ 1 \k=o — is compatible with multiplica- 
tive renormalization. 3 The ultraviolet structure of the 
Faddeev-Popov theory, including asymptotic freedom, is 
preserved. 

These considerations suggest that cutting off the func- 
tional integral at the boundary of the Gribov region f2 
should extend the range of applicability of gauge theory 
toward the infrared, even if ultimately a further refine- 



3 This condition agrees with the Kugo-Ojima confinement criterion 
[2^1 . The relation between the present approach and the Kugo- 
Ojima a ppr oach, heretofore mysterious has been clarified very 
recently l2 r J . 
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ment will be made in which the cut-off is advanced to the 
boundary of the fundamental modular region A. 

Only infrared properties are modified by a cut-off of 
Faddeev-Popov functional integral at the Gribov hori- 
zon, but here the change is dramatic. It implies a linearly 
rising, renormalization-group invariant potential between 
quarks that is transmitted by the bose ghost, which how- 
ever has tensor coupling. Clearly this requires further 
investigation. 

Unresolved open questions are unitarity of the physical 
states, and the relation of confinement to the long-range 
force found here. We remain far from a satisfactory un- 
derstanding the physical states of QCD. 
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APPENDIX A: OTHER EXACT PROPERTIES 

For completeness we note the following relations. The 
mixed V-A propagator is given by 

«(*K(2/)) = -ig(2 7 ) 1 ^{(M- l A x r b (x)Al(y)), 

(Al) 



and the U propagator by 



(U?(x)U* d (y)) = {dxf J ,S bd (M^ 1 ) ac (x, y, A)) 

= 5^5 bd G ac {x-y). (A2) 



where G ac (x — y) is the propagator of the Faddeev-Popov 
ghost. The auxiliary fermi ghosts to and to have the same 
propagator as c-c, and U-U, 



«(aO<(i/)> = S Xfl S bd G ac (x - y). (A3) 



One may also show in the infrared limit 



D AA (k)D AV (k) - D 2 AV {k) = 0(1). (A4) 
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